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Symmetry protected topological (SPT) states are bulk gapped states with gapless edge excitations 
protected by certain symmetries. The SPT phases in free fermion systems, like topological insula- 
tors, can be classified by the K-theory. However, it is not known what SPT phases exist in general 
interacting systems. In this paper, we present a systematic way to construct SPT phases in inter- 
acting bosonic systems, which allows us to identify many new SPT phases, including three bosonic 
versions of topological insulators in three dimension and one in two dimension protected by particle 
number conservation and time reversal symmetry. Just as group theory allows us to construct 230 
crystal structures in 3D, we find that group cohomology theory allows us to construct different 
interacting bosonic SPT phases in any dimensions and for any symmetry groups. In particular, 
we are going to show how topological terms in the path integral description of the system can be 
constructed from nontrivial group cohomology classes, giving rise to exactly soluble Hamiltonians, 
explicit ground state wave functions and symmetry protected gapless edge excitations. 

PACS numbers: 71.27.+a, 02.40.Re 



We used to believe that different phases of matter are 
different because they have different symmetries. Re- 
cently, we see a deep connection between quantum phases 
and quantum entanglement^^ which allows us to go be- 
yond this framework. First it was realized that even 
in systems without any symmetry there can be distinct 
quantum phases - topological phases'"'^ due to different 
patterns of long-range entanglement in the states.^ For 
systems with symmetries, difference in long-range entan- 
glement and in symmetry still lead to distinct phases. 
Moreover, even short-range entangled states with the 
same symmetry can belong to different phases. These 
symmetric short-range entangled states are said to con- 
tain a new kind of order - symmetry protected topolog- 
ical (SPT) order. ^ The SPT phases have symmetry pro- 
tected gapless edge modes despite the bulk gap, which 
clearly indicates the topological nature of this order. On 
the other hand, the gapless edge modes disappear when 
the symmetry of the system is broken, indicating that 
this is a different type of topological order than that 
found in fractional quantum Hall systems^^'^^ whose edge 
modes cannot be removed with any local perturbation.^^ 
Also, SPT orders have no factional statistics or fractional 
charges, while intrinsic topological orders from long range 
entanglement can have them. The discovery of SPT or- 
der hence greatly expands our original understanding of 
possible phases in many-body systems. 

One central issue is to understand what SPT phases 
exist and much progress has been made in this re- 
gard. The first system with SPT order was discovered 
decades ago in spin-1 Haldane chains. The Haldane chain 
with antiferromagnetic interactions was shown to have a 
gapped bulk^^'^^ and degenerate modes at the ends of the 
chain ^^^^ which are protected by spin rotation or time 
reversal symmetry of the system.^' This model has been 
generalized, leading to a complete classification of SPT 



orders in one dimension. Topological insulator 
with time reversal symmetry and particle number con- 
servation protected gapless edge modes provides the first 
example of SPT order in higher dimensions. The non- 
interacting nature of fermions in this system allows a 
complete classification of such kind of SPT orders. ^^'^^ 

However, the understanding of SPT orders in free 
fermion systems is not enough, because particles in many- 
body systems do interact and interaction could dra- 
matically change the phase diagram obtained for free 
fermion systems. Examples have been found in one^^ 
and two^^^^ dimensions where different free fermion SPT 
phases become the same with interactions. Moreover, 
new interacting SPT phases can exist which cannot be 
realized by free fermion systems. Therefore, to have a 
complete understanding of SPT order, we need to an- 
swer the following question: what SPT phases exist in 
general interacting systems, with topological features sta- 
ble under any symmetric interaction as long as no phase 
transition occurs? 

In this paper, we answer this question by presenting 
a systematic construction of SPT phases in interacting 
bosonic systems in any dimension and for any symmetry. 
The stability of the order in the constructed model can 
be proven under any type of symmetric interaction (at 
least in one and two dimension). In one dimension, our 
construction reproduces the classification result already 
known. While no bosonic SPT phases were known pre- 
viously in two and higher dimensions, our construction 
leads to the discovery of many new SPT phases in such 
systems, which are summarized in table I. As listed in 
the first row, we find one kind of bosonic topological in- 
sulator in 2D and three kinds in 3D with boson number 
conservation symmetry U{1) and time reversal symme- 
try Zj. If boson numbers are allowed to fluctuate but 
time reversal symmetry is preserved, then we find one 
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TABLE I: (Color online) SPT phases in d-spatial dimen- 
sions protected by some simple symmetries (represented by 
the symmetry groups). Here Zi means that our construc- 
tion only gives rise to the trivial phase. Zn means that the 
constructed non-trivial SPT phases plus the trivial phase are 
labeled by elements in Zn- Z2 represents time reversal sym- 
metry, U(l) represents boson number conservation symme- 
try, Zn represents cyclic symmetry, etc . The red row is for 
bosonic topological insulators and the blue row is for bosonic 
topological superconductors. 

kind of bosonic topological superconductor in every odd 
spatial dimension, as listed in the second row. More gen- 
erally, for systems in d spatial dimension and with sym- 
metry of group G, we are going to present a way to write 
down a quantized topological term in the path integral 
in d + 1 dimensional space-time based on the nontrivial 
group cohomology of G. From the path integral, we can 
find the ground state wave function, identify the gapless 
edge states and understand how symmetry protects the 
gaplessness of the edge and hence the SPT order even 
against strong interactions. For simplicity, we are go- 
ing to first present the construction in detail for the one 
dimensional Haldane chain, demonstrate the emergence 
of its SPT order from the nontrivial group cohomology 
of time reversal symmetry and then generalize to higher 
dimensions and to all other symmetries. 

The fixed point ground state wave-function of the 
Haldane chain^ takes a simple dimer form as shown in 
Fig.l, where each site contains two spin 1/2 's which are 
connected into singlet pairs | tiii+i) ~ I iiti+i) be- 
tween neighboring sites. Time reversal symmetry acts 
as M(T) = icFyK on each spin 1/2, where K is complex 
conjugation and (jy is the y component of the spin opera- 
tor. The wave-function is invariant under the symmetry 
action. Note that for each spin 1/2 M(T)^ = — / while 
on each site with two spins (M(T) (g) M(T))^ = /. So 
the states on each site form a representation of Zj, the 
symmetry group generated by time reversal symmetry. 
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FIG. 1: (Color online) Dimer form of the ground state wave- 
function in Haldane chain. Each site (big oval) contains two 
spin 1/2's (small dot), which are connected into singlet pairs 
(connected dots) between neighboring sites. 

The wave-function on a closed chain is the gapped 



ground states of anti-ferromagnetic Heisenberg interac- 
tions between each pair of spin 1/2's, H = J^i^I " 
(T^^^ where a\ and <t[ are spin operators for the left and 
right spin 1/2 on each site respectively. Obviously, the 
Hamiltonian is time reversal invariant. Therefore we have 
a situation where the ground state does not break any 
symmetry of the system, yet the system is far from a triv- 
ial phase which becomes evident when we put the system 
on an open chain. When the chain is open, the dangling 
spin 1/2 at each end forms a nontrivial projective repre- 
sentation of Z2 with M{TY = which does not allow 
a one-dimensional representation (for definition of pro- 
jective representations, see appendix A). Therefore, the 
degeneracy of the edge state is robust under any pertur- 
bation as long as time reversal symmetry is preserved. 

Indeed, the nontrivial projective representation on the 
edge is the key to the existence of SPT order in the Hal- 
dane chain and applies in general to all one dimensional 
bosonic SPT phases. ^^'^^ A generic way to write projec- 
tive representations of group G is to define symmetry 
operations M{g){g G G) on group element labeled states 
1^0) (^0 ^ G) as 

M{g)\go) = 4'\go.g-'g''.gn\ggo).g'' = E (i) 

where the 2-cocycle ^^2(^0? ^i? 5^2) is a function from three 
group elements to a U{1) phase factor satisfying 

^l^^\go^gi^g2) = ^2{ggQ,ggi,gg2), g ^G (2) 

. ^2(^1,^2,^3)^2(^0,^1,^3) 

and — r — r 1 (3) 

^2(^0,^2,^3)^2(^0,^1,^2) 

s{g) = 1 if ^ is unitary and 5(^) = — 1 if ^ is antiu- 
nitary. (We prove this in appendix A.) The projective 
representation on the edge of the Haldane chain is of ex- 
actly this form if we relabel the spin states with group 
elements. The time reversal symmetry group contains 
two elements = {E^ T} with T o T = E. For the 
left spin 1/2 on each site, label |t)/|l) as \E)/\T) and 
for the right one, label W)/\i) as \T)/-\E). The total 
wave-function becomes 

\^ain.er) = \{\VTU) + \ElEU,) (4) 

i 

Time reversal symmetry on the edge spin then acts as 
M{T)\E) = -\T) and M{T)\T) = \E), which indeed 
takes the form in Eq.(l) with the 2-cocycle V2 of time 
reversal symmetry given by 

^2(^,r,^) = ^2(r,^,r) = -i (5) 

all other i^2 (^o , ^i , ^2 ) = 1, G 

This essential feature of ID SPT phases can be re- 
produced from a simple construction of the low energy 
effective path integral as we show in the following for the 
Haldane phase. The path integral formulation not only 
provides understanding of the nontrivial physics from a 
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FIG. 2: (Color online) (a) A branched triangularization 
of space-time. (For details of braching see appendix D) 
(b) A tetrahedron - the simplest discrete closed surface. 
Y[i^^^^'^ (gi, gj , Qk) = 1 on a tetrahedron is guaranteed by 
Eq.(3). Note that S123 = sois = 1 and S023 = S012 = —1- (c) 
Discretized space-time manifold Mext on an open disk with 
boundary manifold M. gi G M, ^* is in the interior of Mext- 



bulk perspective but also forms the basis for the gener- 
alization of SPT phases into higher dimensions. 

To write the path integral for disordered phases, we 
first realize that the field of the system fluctuates strongly 
at all length scales and the low energy effective theory has 
no continuous limit. Therefore, we choose to write the 
path integral on discrete space-time. Imagine that we dis- 
cretize the (1 + 1)D space time of the Haldane chain with 
a branched triangularization as shown in Fig. 2 (a). On 
each vertex of the space-time complex, we put a G . 
Time reversal acts as complex conjugation K together 
with a mapping from Qi to Tgi. The path integral is 
then written as 

{9i} 

e-S(te})= -Q 4^^\g^,g^,g^) (6) 

{ijk} 

|G| is the number of elements in G(\G\ = 2 for the time 
reversal symmetry), Ny is the number of vertices in the 
complex. Also, Sijk — ±1 depending on the orienta- 
tion of the triangle. Compared to the continuous for- 
mulation of path integral Z = J Dg q- f dxdr /:[gix,T)] 
gi corresponds to the field g{x^t) and XI {^} corre- 
sponds to the path integral J Dg. Q-^iisn}) the 
action-amplitude on the discretized space-time that cor- 
responds to e~ J' ^ ^^^^^'^''^ of the continuous for- 
mulation and ^2^^ {di^dj^dk) corresponds to the action- 
amplitude Q~ hi,3,k) ^^^^ ^^9{.x,t)] ^ single triangle 
(z, j, /c).The above term can be formally regarded as the 
discrete symmetry group and discrete space time general- 
ization of topological ^-term for an 0(3) non-linear sigma 
model. (For details see appendix C.) 

From the properties of i^2{go^ gi^ g2)^ it can be checked 
that the path integral as defined is symmetric under time 
reversal symmetry with an action amplitude that is in a 
fixed point form and is always equal to 1 on a closed sur- 
face. First because i^2{Tgo,Tgi,Tg2) = i^^^(^o, ^i, ^2) 
(see Eq.(2)), the path integral is invariant under time 
reversal. Secondly, Eq.(3) gives rise to a renormaliza- 



tion flow under which the form of the action amplitude 
remains invariant (for details see appendix D). Finally, 
Eq.(3) also guarantees that the action amplitude on a 
closed space time surface is always 1. This can be easily 
verified on the simplest discrete closed surface - a tetra- 
hedron, as shown in Fig. 2(b). More complicated closed 
surfaces are obtained by putting tetrahedrons together, 
hence the action amplitude will always be 1. Therefore, 
we have constructed a quantized topological term for the 
path integral description of (1 + 1)1^ systems with time 
reversal symmetry on discrete space time. 

But how do we know that this quantized topological 
term describes the SPT order in Haldane chain? This 
can be made explicit by finding the ground state wave- 
function of the system from the quantized topological 
term. The ground state wave-function can be obtained 
by imaginary time evolution from time —00 until time 0. 
In our formulation, this is equivalent to imaginary time 
path integral on a space-time geometry with a boundary 
(at time 0). Denote the boundary as M and the whole 
manifold (a disk) as Mext- As we are considering a fixed 
point path integral, it does not matter how big the inte- 
rior of Mext is and we can reduce it, for example, to just 
one point as shown in Fig. 2 (c). 

To obtain the ground state wave- function, we fix the 
degrees of freedom {gi}M on M and choose an arbitrary 
configuration, for example ^* = for degrees of freedom 
in the interior of Mext- We find 



(7) 



where Yl- is product over all triangles on Mext and for 
simplicity of notation we have chosen all triangles to be 
oriented clockwise. The wave-function on M does not 
depend on the choice of g* . Note that time reversal acts 
as complex conjugate K together with a change of basis 
1^) ~^ I'T) ~^ \^) on each g^ and the wave-function 
is invariant under this action. 
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FIG. 3: (Color online) Duality transformation between wave- 
functions in Eq.(4) and Eq.(7). 



This form of wave-function does not look exactly the 
same as the dimer form shown in Fig.l, but it will be 
after a duality transformation as shown in Fig. 3. First 
expand each g^ into two degrees of freedom and h\_^i 
such that = h\_^i = gi and the amplitude of each 
configuration in the wave-function remains unchanged, 
^{{K = /^i+i = 9i}) = UiM9i^9i+u9*)' Then com- 
bine h\ and into one site and apply a change of basis 
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on each site 

\hi, hi)' = V2{hi K,g*)\hi hi) = M9i-u9i, 9*)\hl, K) 

The amphtude of ah configurations in the new basis be- 
comes 1, ^'{{Ki = h\_^i = gi}) = 1 which can be equiva- 
lently written as a product of dimers between neighboring 
sites = Yl- Y^g, \h\ = Qi, h\^^ = gi). In this way we 
have mapped each degree of freedom gi into a dimer and 
the total wave- function takes the same form as Eq.(4). 

It is also important to check that time reversal sym- 
metry acts as in Eq.(l) on the edge degree of freedom in 
the ground state. To see this, note that after the duality 
transformation, time reversal symmetry acts separately 
on each site and is composed of complex conjugation K 
together with unitary transformation on the basis states 



MU9)\h\,K)' 



vr\h\,hl,g*) 
V2{gh\,ghl,g*) 

4'\hi9-'9%9*), u 



\9higKy 
\9K,9K)' 



which factorizes on the two degrees of freedom and acts 
in the same way as in Eq.(l) on the edge. Therefore, 
our quantized topological term Eq.(6) provides a proper 
description of the SPT order in the Haldane chain. 

The advantage of reformulating the Haldane chain in 
the path integral language is that it can be easily general- 
ized to all spatial dimensions d and all symmetry groups 
G. The two cocycles ^^2(^0, ^i, ^2) used in the construc- 
tion have higher dimensional analogues - the {d ^ 1)- 
cocycles ^^(^+1(^0, ^d+i), which are maps from d + 2 
group elements to a U{1) phase factor and satisfy 

^d+i(^o,^i,...,^d+i) = Vd+i{9gQ,ggi,"',g9d+i),g ^ G 

(^0,..,^i-l,^i+l,...,^d+2) = 1 

i=0 

A quantized topological term in the path integral of a 
{d + 1)D system can be written in terms of the {d + 1)- 
cocycles by (1) discretizing the {d-\-l)D space time with 
triangularization (triangle in {l-\-l)D and tetrahedron in 
(2-h 1)D etc ); (2) putting group element labeled degrees 
of freedom onto the vertices (3) assigning action ampli- 
tude to each simplex with the corresponding cocycle. The 
path integral then takes the form 



\G\ 



n ^7+T"i9i,9j,-,9k) (8) 

{gi} {ij...k} 



where = ±1 depending on the orientation of the 

simplex ij...k. Similar to the (1 + 1)D case, it can be 
shown that the path integral is symmetric under sym- 
metries in group G, the action amplitude e"*^*^^^^^^ = 
ll{ij...k} ^d+i^idiigji "',gk) is in a fixed point form and 
is quantized to 1 on a closed manifold. 



The ground state wave-function can be obtained from 
the action amplitude on a open geometry as discussed be- 
fore ^{{gi}M) n{i...j*}^2(^i,...,^j,^*), where {^Jm 
is on M and ^* is inside Mext- j*} is the prod- 

uct over all simplices. An exactly soluble Hamilto- 
nian can be constructed to realize this state as the 
gapped ground state. To see this, note that ^ can be 
mapped through a local unitary transformation U = 

Y{{i...j^]^2^{gi^---^gj^g^)\gi---gj)\gi---gj\ a total prod- 
uct state ^o({^i}M) = 1- As ^0 is the gapped 
ground state of i^o = - E/c \^k){(t>kl Vt>k) = T^g.eG 
^ is the gapped ground state of H = WHoU = 
—U^(l)k){4^k\U, which is local and has symmetry G. 
The nontrivial SPT order of the system can be seen 
explicitly from its boundary. The path integral of degrees 
of freedom on the boundary can be obtained by putting 
the topological term on an open geometry as shown in 
Fig. 2 (c) for {1^1)D. Note that the manifold M now 
corresponds to the space time manifold of the boundary 
degrees of freedom. The path integral for the boundary 
then reads 



-N 



E n ^d+ri9i,-,9j,9*) (9) 

{9i}M {^---J*} 



where {gi}M is the field on M and g* is inside Mext- The 
action amplitude for each {gi}M configuration does not 
depend on the choice of 

This term can be thought of as a discretized version of 
the Wess-Zumino-Witten(WZW) term^^'^^ in nonlinear 
<j-models because 1. it is a path integral of {d — 1) -\- 1 
dimensional systems written on an extended {d + 1)D 
manifold with a boundary 2. the action amplitude does 
not depend on how the extended field in the interior of 
the {d -\- 1)D manifold is chosen 3. its field takes value 
in a group G and the path integral is invariant under the 
action of ^ G G. On the other hand, this term is more 
general than the original continuous WZW term because 
it applies to discrete groups like Z2 while the contin- 
uous WZW term only works for continuous groups. We 
expect that the boundary theory described by such a dis- 
cretized WZW term will be gapless as long as symmetry 
is not broken, similar to systems described by continuous 
WZW terms. Actually this has been firmly established in 
(1 + 1)D and (2 + 1)D. In (1 + 1)D, as we have discussed 
with the example of time reversal symmetry, symmetry 
action on the edge degree of freedom does not have a one- 
dimensional representation, therefore the edge state will 
always be degenerate. In (2 + 1)D^ we have proved rigor- 
ously using the tool of matrix product unitary operator 
that the symmetry action on the boundary requires it to 
be gapless as long as symmetry is not broken. ^^'^^ There- 
fore, the boundary of the systems we constructed carries 
gapless excitations protected by certain symmetry, which 
reflects the nontrivial SPT order of the system. 

In conclusion, we have presented a systematic con- 
struction of symmetry protected topological phases in 
bosonic systems in any spatial dimension d and with any 
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symmetry of group G. In particular, we wrote down a 
quantized topological term in the path integral of the 
system based on nontrivial {d + l)-cocycles of the sym- 
metry group. From the quantized topological term we 
found that the ground state wave-function is short range 
entangled and the boundary of the systems carries gap- 
less excitations protected by the symmetry of the system, 
which clearly demonstrates the nontrivial SPT order of 
the system. The number of nontrivial cocycles, hence 
the number of nontrivial SPT phases, for some simple 
symmetry groups are summarized in table I. Our con- 
struction is non-perturbative and works for strongly in- 
teracting bosonic systems. 

This research is supported by NSF Grant No. DMR- 
1005541 and NSFC 11074140. Z.C.G. is supported by 
NSF Grant No. PHY05-51164. 



Appendix A: Projective Representation 

Matrices u{g) form a projective representation of sym- 
metry group G if 

u{gi)u{g2) = uj{gi,g2)u{gig2), gi,g2 ^ G. (Al) 

Here uj{gi^g2) ^ ^(1) and uo{gi^g2) 7^ 1, which is called 
the factor system of the projective representation. The 
factor system satisfies 



{g2, 93)^(91^ 9293) =uj{gi,g2)^{9i92,93)^ (A2) 



for all ^1,^2,^3 ^ G, where s{gi) = 1 if ^1 is unitary and 
s{gi) = -1 if gi is anti-unitary. If cj(^i, ^2) = 1, V^i,^2, 
this reduces to the usual linear representation of G. 

A different choice of pre-factor for the representation 
matrices u\g) = /3{g)u{g) will lead to a different factor 
system uj\gi,g2): 



^'{91^92) 



P{gi)f3<^^\g2) 
P{9i92) 



uj{gi,g2). (A3) 



We regard u'{g) and u{g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems oj'{gi^g2) and uo{gi^g2) as be- 
longing to the same class uo. 

Suppose that we have one projective representation 
Ui{g) with factor system (^1(^1,^2) of class uoi and an- 
other U2{g) with factor system 002(91^92) of class CJ2, ob- 
viously Ui{g)^U2{g) is a projective presentation with fac- 
tor system ^2)^2(^1, ^2)- The corresponding class 
uj can be written as a sum uji + co'2. Under such an addi- 
tion rule, the equivalence classes of factor systems form 
an Abelian group, which is called the second cohomology 
group of G and is denoted as 1-C^[G^ The identity 
element 1 G H^[G, 11(1)] is the class that corresponds to 
the linear representation of the group. 

^{91^92) can be equivalent ly expressed as 1^2(90^ 9i^ 92) 
with correspondence 



where E is the identity element in the symmetry group 
G. 1^2(90,91^92) satisfies 



sig) 



(90,91,92) = ^2(990,991,992), 9 (A5) 



and the condition Eq.(A2) satisfied by oj(gi^g2) becomes 
^2(91, 92, 93)^2(90, 91, 93) 



^2(90, 92, 93)^2(90, 91, 92) 



1 



(A6) 



Similarly, P(gi) can be equivalently written as vi(E^gi) 
where 



^l^^\9o,9i) = ^1(990,991) 



(A7) 



^(91,92) = V2(E,gi,gig2) 



(A4) 



The equivalence relation Eq.(A3) between cj(^i,^2)'s 
then becomes 

// A ^1(90,91)^1(91,92) ( \ f\Q\ 

1^2(90,91,92) = 7 ^ v(gQ,gi,g2). (A8) 

In the following we will show M(g) defined in Eq.(l) 
forms a projective representation of symmetry group G. 
It is easy to check that: 

M(gi)M(g2)\gQ) 
=M(g,)4'^\g,,g^'g\g^)\g2g,) 

=4''\9290,9i'9\9l4''''\90,92 '9\9l\9i9290^ 

(A9) 

with g* = E. We note that in the last line, the sign 
factor s(gig2) arises because the complex conjugate K 

could have a nontrivial action on i'2^^'^\9o, 92^ 9* , 9*) if 
gi is an antiunitary. According to the 2-cocycle condition 
Eq. (A6), we have: 

M(gi)M(g2)\go) 
=^2 (9* ,9i9*, 9i929*)^2 (919290 , 9* , 91929* ) \9i929o) 

=^2(9"^ , 919"" , 9i929l^2^^'^^\9o, (9i92)~^9*,9*)\9i929o) 
=U2(E, gi, gig2)M(gig2)\go) 

=uJ2(9i,92)M(gig2)\go) (AlO) 

Now it is clear that M(g) defined in Eq.(l) forms a pro- 
jective representation of the symmetry group G. How- 
ever, M(g) is usually reducible. It can be reduced to 
a direct sum of several irreducible projective representa- 
tions which belong to the same class. Only irreducible 
projective representations describe the edge states of gen- 
eral ID SPT phases. 



Appendix B: Group cohomology 

The above discussion on the factor system of a projec- 
tive representation can be generalized which gives rise to 
a cohomology theory of groups. In this section, we will 
briefiy describe the group cohomology theory. 
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For a group G, let M be a G-module, which is an 
abehan group (with multiphcation operation) on which G 
acts compatibly with the multiplication operation {ie the 
abelian group structure) on M: 

g • {ah) = {g • a){g -b), g e a,b e M. (Bl) 

For the cases studied in this paper, M is simply the U{1) 
group and a a /7(1) phase. The multiplication operation 
ah is the usual multiplication of the U{1) phases. The 
group action is trivial g-a = a{geG,ae U{1)) if g is 
unitary and ^ • a = a* if ^ is anti-unitary. 

Let z^n(^o, "'i9n) be a function of (n+1) group elements 
whose value is in the G-module U{1). In other words, 
Un : G^^+i) ^ U{1). Vn satisfies 

<^^H^o,^i,...,^n) = Vn{99Q,99i,'",99n),9 e G (B2) 

We will call such a map an n-cochain: Let 
C^{G, U{1)) = {i^n} be the space of all n-cochains. Note 
that C^(G, U{1)) is an Abelian group under the function 
multiplication v'^{gQ, gn) = i^n(^o, ^nX(^o, ^n). 
We define a map from C^[G, U{1)] to C^+i[G, U{1)]: 

{dn^n){9Q,9li "'i9n+l) 
n+1 

= n ^n"^^'(^0,..,^i-l,^i+l,...,^n+l) (B3) 

Let 

B"(G, U{1)) = {VnK = dn-lI^n-U ^ C^-\G, U{1))} 

and 

Z-(G, U{1)) = {Vn\dnyn = 1, ^ C^(G, U {!))} 

S^(G,/7(1)) and 2:^(G,/7(1)) are also Abelian groups 
which satisfy S^(G,/7(1)) C ^^(G, [/(!)) where 
Si(G, U{1)) = {!}. The n-cocycle of G is defined as 

H^(G, /7(1)) = Z^(G, /7(1))/S^(G, U{1)) (B4) 

Let us discuss some simple cases. From 

{diiyi){go,gi,g2) = i^iigo, gi)iyi{gi, g2) /M9o^ 92) 
we see that 

Z^(G, /7(1)) = {ui\ui{go,gi)iyi{gi,g2) = ^^1(^0,^2)}. 

Since B\G,U{1)) = {1} is trivial, n\G,U{l)) = 
Z^{G^U{1)). If we define a{g) = vi{E^g)^ it is easy to 
see that the first-cocycle is related to the one-dimensional 
representations of the group G. 
From 

(g^2Z^2)(^0,^1,^2,^3) 

= 1^2(91^ 92, 9s)M9o, 91, 93) /M90, 92, 93)^90,91,92) 



we see that 
Z\GMl)) = {^2\ 

^^2(91, 92, 93)^^2(90, 91,93) 1^2(90, 92, 93)^^2(90, 91, 92)}' 
and 

B\G,U{1)) = {U2\ 

^^2(90,91,92) = 1^1(90, gi)M9i, 92) /j^i{9o, 92)}' 

The 2-cocycles ^^(G, [/(!)) = {G , U {!)) / {G , U {!)) 
classify the projective representations discussed in sec- 
tion A. 
From 

{d3J^3){9o,9i,92,93,9^) 

^3(91, 92, 93, 94)^^3(90, 9i, 93, 94)^3(90, 9i, 92, 93) 



^3{9o,92,93,9a)^3{9o,9i,92,9a) 



(B5) 



we see that 



Z\G,U{l)) = {v,\ 



'^3(50,52,53,54)1^3(50,51,52,54) 



= !}■ 



and 



B\G,U{l)) = {v,\ 

( N J^2{9i, 92, 93)^90,91,93). 

m9o,9i,92,93) = rj 

^^2(^0, ^2, ^3)^2(^0, ^1,^2) 

which give us the 3-cocycle H^(G, [/(I)) = 
Z3(G,[/(1))/S3(G, [/(!)). 



Appendix C: Relationship with the topological-^ 
term of the 0(3) non-linear sigma model 

In the section, we will show the amplitudes Eq.(6) can 
be formally regarded as the discrete group and discrete 
space time generalization of topological term of non- 
linear sigma model. We will first review that the Haldane 
phase (a non-trivial ID SPT phase) is described by a 27r- 
quantized topological term in continuous non-linear a- 
model.^^ However, such kind of 27r-quantized topological 
terms cannot describe more general ID SPT phases. We 
argue that to describe SPT phases correctly, we must 
generalize the 27r-quantized topological terms to discrete 
space-time. 

Before considering a spin-1 chain, let us first consider a 
(0 + 1)D non-linear <j-model that describes a single spin, 
whose imaginary-time action is given by 

■ is [ dtd^ n(t,0 • [dtn{t,0 X d^n{t,^)] 

(CI) 
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where n(t) is an unit 3d vector and we have assumed 
that the time direction form a circle. The second term 
is the Wess-Zumino-Witten (WZW) term.^'''^^ We note 
that the WZW term cannot be calculated from the field 
n(t) on the time-circle. We have to extend n{t) to a 
disk D'^ bounded by the time-circle: n{t) (see 
Fig. 4). Then the WZW term can be calculated from 
n(t, ^). When 2s is an integer, WZW terms from different 
extensions only differ by a multiple of 2i7r. So e"*^ is 
determined by n{t) and is independent of how we extend 
n{t) to the disk D'^. 




FIG. 4: If we extend n(t) that traces out a loop to n(t, ^) that 
covers the shaded disk, then the WZW term f^2 dtd^ n(t, ^) • 
[dtn{t,^) X d^n{t,^)] corresponds to the area of the disk. 

The ground states of the above non-linear a-model 
have 25 + 1 fold degeneracy, which form the spin- 5 rep- 
resentation of 5*0 (3). The energy gap above the ground 
state approaches to infinite as ^ ^ oo. Thus a pure 
WZW term describes a pure spin-5 spin. 




FIG. 5: (a) The topological term W describes the number 
of times that n{x,t) wraps around the sphere (as we change 
t). (b) On an open chain x G [0,L], the topological term W 
in the (1+1)D bulk becomes the WZW term for the end spin 
TiLit) = n(L,t) (where the end spin at x = is hold fixed). 

To obtain the action for the SO {3) symmetric antifer- 
romagnetic spin-1 chain, we can assume that the spins 
Si are described by a smooth unit vector field n{x^t): 
Si = {—yniia^t) (see Fig. 5b). Putting the above single- 
spin action for different spins together, we obtain the fol- 
lowing (1+1)D non-linear a-model"^^ 

S = I dxdt — {dn{x,t)f ^lOW, (9 = 27r, (C2) 

where W = (47r)~-^ / dtdx n{t, x) • [dtn{t, x) x dxn{t, x)] 
and iOW is the topological term.^^ If the space-time 
manifold has no boundary, then e~^^^ = 1 when = 



mod 27T. We will call such a topological term - a 
27r-quantized topological term. The above non-linear a- 
model describes a gapped phase with short range corre- 
lation and the SO {3) symmetry, which is the Haldane 
phase. In the low energy limit, g flows to infinity and 
the fixed-point action contains only the 27r-quantized 
topological term. Such a non-linear <j-model will be 
called topological non-linear a-model. 

It appears that the 27r-quantized topological term has 
no contribution to the path integral and can be dropped. 
In fact, the 27r-quantized topological term has physical 
effects and cannot be dropped. On an open chain, the 
27r-quantized topological term 27riW becomes a WZW 
term for the boundary spin nL{t) = n{x — L^t) (see 
Fig. 5).^^ The motion of is described by Eq.(Cl) 
with s = 1/2. So the Haldane phase of spin-1 chain has 
a spin-1/2 boundary spin at each chain end!^^'^^'^^ 

We see that the Haldane phase is described by a fixed 
point action which is a topological non-linear a-model 
containing only the 27r-quantized topological term. The 
non-trivialness of the Haldane phase is encoded in the 
non trivially quantized topological term.^^'^^ From this 
example, one might guess that various SPT phases can 
be classified by various topological non-linear a- models, 
and thus by various 27r-quantized topological terms. But 
such a guess is not correct. 

This is because the fixed-point action (the topological 
non-linear a-model) describes a short-range-correlated 
state. Since the renormalized cut-off length scale of the 
fixed-point action is always larger than the correlation 
length, the field n{x^t) fluctuates strongly even at the 
cut-off length scale. Thus, the fixed-point action has 
no continuum limit, and must be defined on discrete 
space-time. On the other hand, in our fixed-point action, 
n(x, t) is assumed to be a continuous field in space-time. 
The very existence of the continuum 27r-quantized topo- 
logical term depends on the non-trivial mapping classes 
from the continuous space-time manifold to the con- 
tinuous target space S'^. It is not self consistent to use 
such a continuum topological term to describe the fixed- 
point action for the Haldane phase. 

As a result, the continuum 27r-quantized topological 
terms fail to classify bosonic SPT phases. For exam- 
ple, different possible continuum 27r-quantized topologi- 
cal terms in Eq.(C2) are labeled by integers, while the 
integer spin chain has only two gapped phases protected 
by spin rotation symmetry: all even-integer topological 
terms give rise to the trivial phase and all odd-integer 
topological terms give rise to the Haldane phase. Also 
non-trivial SPT phases may exist even when there is no 
continuum 27r-quantized topological terms (such as when 
the symmetry G is discrete). 

However, the general idea of using fixed-point actions 
to classify SPT phases is still correct. But, to use 27r- 
quantized topological terms to describe bosonic SPT 
phases, we need to generalize them to discrete space- 
time. In the following, we will show that this indeed can 
be done, using the (1+1)D model (C2) as an example. 
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A discrete (1+1)D space-time is given by a branched 
triangularization (see Fig. 2 (a)). Since S = J dxdtL, 
on triangularized space-time, we can rewrite 



(jii.nj.rik) = e 



— r. dxdt L 



e/7(i), 



(C3) 



where dxdt L is the action on a single triangle. We see 
that, on discrete space-time, the action and the path inte- 
gral are described by a 3- variable function i^(ni, n^, n/c), 
which is called action amplitude. The 5*0 (3) symmetry 
requires that 

v{gni,gnj,gnk) v{ni,nj,nk), g e S0{3). (C4) 

In order to use the action amplitude i^^^*'-^'^^(n^, rij, n/e) 
to describe a 27r-quantized topological term, we must 
have ]^ z/^'^*'-^'^^(ni, nj, n/c) = 1 on any sphere. This can 
be satisfied iff Yl i^**^*'-^'^^(ni, n^, n/e) = 1 on a tetrahedron 
- the simplest discrete sphere (See Fig. 2(b)): 



i/(ni, 712, ns)u{no, ni, 713) 
z/(no, 712, 713)^/(710, ni, 712) 



1. 



(C5) 



(Another way to define topological term on discretized 
space-time can be found in Ref. 44.) 

A i/(7io, 7ii, 712) that satisfies Eq.(C4) and Eq.(C5) is 
called a 2-cocycle. If ^^(710, 7ii, 712) is a 2-cocycle, then 



i/'(7lo,ni,7l2) i/(7lo,7li,7l2) 



/i(7li,7l2)/i(no,7li) 
/i(7lo,7l2) 



(C6) 



is also a 2-cocycle, for any /i(7io,7ii) satisfying 
jii{gno,gni) = /i(7io,ni), g e S0{3). Since u{no,n2,ns) 
and ly' {no^n2^ns) can continuously deform into each 
other, they correspond to the same kind of 27r-quantized 
topological term. So we say that i/(7io, 7i2, 713) and 
u^{no^n2^ns) are equivalent. The equivalent classes of 
the 2-cocycles i/(7io, 7i2, 713) give us /7(1)] - the 

2-cohomology group of sphere S'^ with U{1) coefficient. 
H^(6'^, U{1)) classifies the 27r-quantized topological term 
in discrete space-time and with S*^ as the target space. 

Does n^[S^, U{1)] classify the SPT phases with S0{3) 
symmetry? The answer is no. We know that S'^ is just 
one of the symmetric spaces of SO (3). To classify the 
SPT phases, we need to replace the target space S'^ by 
the maximal symmetric space, which is the group it- 
self SO (3) (see Ref. 45 for more discussions). So we 
need to consider discrete non-linear a-model described 
by v{gi,gj,gk), gi.gj.gk e S0{3). Now the 2-cocycle 
conditions becomes 



^{g9i^99j^99k) = iy{9i,gj,gk) ^ U{1), 
iy{gi,g2, 93)^^(90^ 91^93) ^ 
^^(^0, ^2, ^3)^(^0, ^1,^2) 



(C7) 



which defines a "group cohomology" H^[6'0(3), U{1)]. It 
classifies the 27r-quantized topological term for the max- 
imal symmetric space. It also classifies the SPT phases 
with 50(3) symmetry in (1+1)D. 



Appendix D: Branched triangulation and topological 
invariant amplitudes 

As we have shown, the bosonic symmetry protected 
order in d + 1 D can be described by the amplitude: 



~\Q\N^ i 1 ^d-\-l \9i ^9j-,' 



.9k). (Dl) 



where gi^gj^ - - • ^g^ G G are group elements of the sym- 
metry group G and \G\ is the order of G. The d + 1- 
cocycle {gi^ gj^ ■ • • ,^/c) is defined on the branched 

d-\- 1-simplex, with s^j.../e = ±1 uniquely determined by 
the orientation of the corresponding d + 1-simplex. In 
the following, we will show they are topological invari- 
ant amplitudes. Hence, they represent a class of fixed 
point amplitudes of (symmetry protected) topologically 
ordered phases. 

First, we need to give a branching structure to the 
discretized space- time. A branching is a choice of an ori- 
entation of each edge of an n-simplex such that there is 
no oriented loop on any triangle. For example. Fig. 6 (a) 
is a branched 2-simplex and (c) is a branched 3-simplex. 
However, (b) is not an allowed branching because all its 
three edges contain the same orientations and thus form 
an oriented loop, (d) is also not allowed because one of its 
triangle contains an oriented loop. It is easy to check that 
any consistent branched structure can induce a canoni- 
cal ordering for the vertices of n-simplex. Indeed, the 
branching structure will also induce an canonical orien- 
tation of the n-simplex, see Fig. 7. 







(c) 



(d) 



FIG. 6: Examples of allowed ( (a), (c)) and unallowed ((b),(d)) 
branching for 2-simplex and 3-simplex. 

To understand the geometric meaning of the topologi- 
cal invar iance of the above amplitude, let us first discuss 
(1+1)D fixed-point action- amplitude with a symmetry 
group G. For a (1 + 1)D system on a complex with a 
branching structure, a fixed-point action-amplitude (see 
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mmor 





(a) 



FIG. 7: (Color online) (a): A branching structure of a n- 
simplex will induce a canonical ordering for its vertices. For 
example, if a 3-simplex contains a vertex with no incoming 
edge, then we can label this vertex as and canonically la- 
bel the vertices of the rest 2-simplex as v^^v^^v^. Such a 
scheme can be applied for arbitrary n-simplex if n-simplex 
has a canonical label, (b): A branching structure will also in- 
duce a canonical orientation for the n-simplex. For example, 
we can use the right hand rule to determine the orientations of 
surfaces of the tetrahedron, then the orientation of the tetra- 
hedron can be uniquely determined by the surface opposite 
to g° 




FIG. 8: (Color online) The graphic representation of the 
action- amplitude e"'^*^"'^^*^'* on a complex with a branching 
structure represented by the arrows on the edge. The vertices 
of the complex are labeled by i. Note that the arrows never 
for a loop on any triangle. 



Fig. 8) 



-Siigi}) 



{ijk} 



^2{gi, 90,9^)^2^ {9b, 90, 9a) (D2) 

where each triangle contribute to a phase factor 
T^2'^{9ii9ji9k), Y{{ijk} multiply over all the triangles in 
the complex Fig. 8. Note that the first variable gi in 
^2{9ii 9j^ 9k) correspond to the vertex with two out going 
edges, the second variable gj to the vertex with one out 
going edges, and the third variable g^ to the vertex with 
no out going edges. Sijk = ±1 depending on the orienta- 



tion of i ^ j ^ /c to be anti-clock-wise or clock-wise. 

In order for the action- amplitude to represent a quan- 
tized topological ^-term, we must choose 1^2(91^ 9j^ 9k) 
such that 



-Si{9^}) 



H 4'^\g,,g,,gu) = l (D3) 



{ijk} 



on closed space-time complex without boundary, in par- 
ticular, on a tetrahedron with four triangles (see Fig. 
2(b)): 



-S{{9^}) 



n ^2'\9i^9j^9k) 

{ijk] 

^2(^1,^2,^3)^2(^0,^1,^3) 
^2(^0,^1,^2)^2(^0,^2,^3) 



= 1 



(D4) 



Also, in order for our system to have the symmetry gen- 
erated by the group G, its action-amplitude must satisfy 

^-s{{g.}) ^ ^-Si{gg^})^ ^ contains no T 
(^Q-S{{g^})^'^ ^ Q-Si{gg.})^ -f ^ contains one T (D5) 

where T is the time-reversal transformation. This re- 
quires 



^'2^^\9i^9j^9k) = M99i^99j^99k)' 



(D6) 



Eq.(D5) and Eq.(D6) happen to be the conditions of 2- 
cocycles 1^2(90^ 9i^ 92) of G. Thus the action-amplitude 
Eq.(D2) constructed from a 2-cocycle ^^2(^0, ^1, ^2) is a 
quantized topological ^-term. 

If ^2(90^ 92^ 93) satisfy Eq.(D5) and Eq.(D6), then 



^2(^0,^2,^3) = ^2(^0,^2,^3 



^ Mi(^i, ^2)^1(^0,^1) 
Ml (^0,^2) 



(D7) 



also satisfy Eq.(D5) and Eq.(D6), for any /ii(^o,^i) sat- 
isfying iii{ggQ,ggi) = /ii(^o,^i), ^ ^ G. So 1^2(90^92,93) 
also gives rise to a quantized topological 6>-term. As 
we continuously deform /ii(^o,^i), the two quantized 
topological ^-terms can be smoothly connected. So we 
say that the two quantized topological 6>-terms obtained 
from 1^2(90,92,93) and ^^2(^0, ^2, ^3) are equivalent. We 
note that 1^2(90,92,93) and 1^2(90,92,93) differ by a 2- 
coboundary ^^^3i^g'2)^j^i{go,gi) ^ ^.j^^ equivalence classes 

of 1^2(90,92,93) is nothing but the cohomology group 
H^[G, t/T(l)]- Therefore, the quantized topological 9- 
terms are classified by H^[G, /7t(1)]- 

We can also show that Eq.(D2) is a fixed-point action- 
amplitude from the cocycle conditions on V2(9i, 9j, 9k)' 
The cocycle condition Eq.(D4) induces two renormaliza- 
tion moves in the discretized manifold. Fig. 9 represents 
a 2 ^ 2 moves: 

1^2(90,91,93)1^2(91,92,93) = 1^2(90,91,92)1^2(90,92,93) 

(D8) 
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FIG. 9: (Color online) Graphic representation of 

1^2(90,91,92)1^2(90,92,93) = 1^2(91, 92, 93)1^2(90, 91, 93) The ar- 
rows on the edges represent the branching structure. 




FIG. 10: (Color online) Graphic representation of 

1^2(91,92,93) = 1^2(90, 91, 92)1^2(90, 92, 93)1^2 ^ (90, 91, 93)- The 
arrows on the edges represent the branching structure. 

and Fig. 10 represents a 1 ^ 3 move: 

^2(^1,^2,^3) ^2(9^,91,92)^2(90,92,9^)^2^(9^,91,9?) 

(D9) 

By using these two moves, different triangularization of 
the space-time can be mapped into each other without 
changing the action amlitude Eq.(Dl). Therefore, the 
action-amplitude is a fixed point form. 

Geometrically, these two moves can be obtained by 
projecting a 3 simplex onto two two-dimensional surfaces. 



Fig. 11 (a) gives the 2^2 move and Fig. 11 (b) gives the 
1^3 move. Note that the projection from opposite 
directions will induce opposite chiralities for the bound- 
ary of the tetrahedron, that's why we need to change the 




(a) (b) 



FIG. 11: (d2T^2)(9o, 91, 92, 93) can be represented as the 
boundary of a 3-simplex. (a) and (b) correspond to two dif- 
ferent basic moves of 2-simplexes. 

chiralities of the triangular in one side of basic moves. 
Such a change correspond to inverse V2 when we move 
it from left side to right side of the above equation, 
which is consistent with our rules of algebra. Similar 
process also works in higher dimensions. Different moves 
in d-dimension can be obtained from projecting a d + 1- 
simplex in different ways. Such moves induces a renor- 
malization flow in the discrete space-time under which 
the action- amplitude Eq.(Dl) is a fixed point. 
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